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Abstract 

We study the phase transition from two-flavor nuclear matter to quark matter. A mean field 
model, constructed at the quark level, is used to give the equation of state for nuclear matter, while 
the equation of state for color superconducting quark matter is calculated within the NJL model. 
It is found that at low temperature, the phase transition from nuclear to color superconducting 
quark matter will take place when the density is of order 2.5po - &Po- At zero density, the quark 
phase will appear when the temperature is larger than about 148 MeV. Within the mean field 
treatment, the phase transition from nuclear matter to quark matter is always first order, whereas 
the transition between color superconducting quark matter and normal quark matter is second 
order. 

PACS numbers: 21.65.+f; 25.75.Nq; 12.38.Mh; 12.39.-x 



1 



I. INTRODUCTION 



Hadronic matter is expected to undergo a phase transition to quark matter at high 
temperature and/or high baryon density. Quark matter may exist in the core of neutron star 
where the density is high or it may be produced in the laboratory in heavy ion collisions. 
In quark matter, the presence of a weak attraction between the quarks will result in the 
formation of a condensate of quark pairs. Pairs of quarks cannot be color singlets, and in 
QCD with two flavors of massless quarks, the Cooper pairs favor the formation of a color 3, 
flavor singlet condensate. The color superconductivity of quark matter has been discussed 
in many places in the literature (see for example [l| - 0] arid references therein). 

The phase diagram for chiral symmetry restoration, including possible color supercon- 
ductivity, for quark matter containing just u and d quarks, was discussed in Refs. jE], 0], 
where a rich phase structure was obtained. It was shown that at high temperature, chiral 
symmetry was restored via a second-order phase transition. At low temperature, the phase 
transition was first order and at this time, a new phase with a condensate of Cooper pairs 
appeared. However, we should notice that this phase transition is not a transition from 
hadronic matter to quark matter but a transition from chiral symmetry breaking quark 
matter to chiral symmetry restored quark matter. We know that in the real world, there 
is no quark matter at low density and temperature, rather the chiral symmetry breaking 
quark matter should be replaced by hadronic matter. Therefore, when studying the phase 
diagram, the equations of state (EOS) of hadronic matter are needed. 

To study the properties of hadronic matter we need phenomenological models, since QCD 
cannot yet be used directly. The symmetries of QCD can be used to constrain the hadronic 
interactions and models based on SU(2)l x SU(2)r symmetry and scale invariance have 
been proposed. These effective models have been widely used to investigate nuclear matter 
and finite nuclei, both at zero and at finite temperature @]-[§]. Papazoglou et al. extended 
the chiral effective models to SU(3)l x SU(3)r, including the baryon octet [Io|, As well 



as models based on hadronic degrees of freedom, there are others based on quark degrees of 



freedom, such as the quark meson coupling model [12], [13[], the cloudy bag model [141 ]. the 
Nambu-Jona-Lasinio (NJL) model [l5[ and the quark mean field model [3], etc.. Several 
years ago, a chiral SU (3) quark mean field model based on quark degrees of freedom was 
proposed by Wang et al. [13, lq - In this model, quarks are confined in the baryons by 
an effective potential. The quark-meson interaction and meson self-interaction are based 
on 577(3) chiral symmetry. Through the mechanism of spontaneous symmetry breaking the 
resulting constituent quarks and mesons (except for the pseudoscalars) obtain masses. The 
introduction of an explicit symmetry breaking term in the meson self-interaction generates 
the masses of the pseudoscalar mesons, which satisfy the relevant PCAC relations. The 
explicit symmetry breaking term in the quark-meson interaction gives reasonable hyperon 
potentials in hadronic matter. This chiral SU(3) quark mean field model has been applied 
to investigate nuclear matter [19], strange hadronic matter [17], finite nuclei, hypernuclei 
[Hi], and quark matter [20]. Recently, we improved the chiral SU(3) quark mean field model 
by using a microscopically justified treatment of the center of mass correction [2l[. This new 
treatment has been applied to the study of the li quid -gas phase transition of asymmetric 
nuclear matter as well as strange hadronic matter [22l . 23[. By and large the results are in 
reasonable agreement with existing experimental data. 

As we know, the phase transition between hadronic matter and quark matter can only 
occur at high density or temperature. As the density grows, one expects that the nucleons 
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will begin to overlap. This is expected to be a source of short distance repulsion and 
correlations. Such effects are also encountered in atomic physics and a method of estimating 
the effect was introduced in Ref . (24[ • This was later applied in the context of dense nuclear 
matter by [25[ and we also use this method to estimate the importance of the effect. As we 
will show later, this procedure generates a repulsive interaction between the nucleons and 
increases their chemical potentials in nuclear matter. 

In this paper, we study the phase transition from nuclear matter to quark matter. The 
effect of nucleon finite size is considered in nuclear matter and color superconductivity is 
discussed in quark matter. Since this is the first time we study the volume effect on the phase 
transition in our model, we restrict our investigation to just two flavors, without including 
strange quarks or hyperons. The paper is organized in the following way. Section II sets out 
the formalism of the NJL model which is applied to describe quark matter. In section III 
we introduce the SU (3) quark mean field model which is used to describe nuclear matter 
- although in this work only the non-strange sector of the model is required. The effect of 
finite nucleon size is discussed in this section. The numerical results are presented in section 
IV. Section V is the summary. 



II. NJL MODEL 



In order to study the color superconductivity of quark matter, we choose the NJL model, 
within which it is simple to incorporate the diquark condensate in the Lagrangian. The 
Lagrangian density for the SU (2) NJL model [26] can be written as 



£ x = C {yfa„ - M) < + Gi fcWfyrfl - (r G )J (r a )[ 



-B\ 



where M is the current mass matrix of the u and d quarks, r a (a = 1 — 3) are the Pauli 
matrices. B' is the parameter which will give a reasonable bag constant, B, in quark matter. 
In order to describe the diquark condensate, the interaction of the color 3 fermion bilinears 
should be added. This four-fermion interaction can be written as 



£•>. — Gy 



(2) 



where C denotes the charge conjugation matrix. In the above equations, G\ and G2 are the 
coupling constants. 

From the above Lagrangian, the thermodynamic potential Q, which determines the field 
equations can be obtained. It is convenience to introduce bosonic collective fields into the 
functional integral to make the four fermion interaction into a fermion-boson coupling. We 
follow Ref. 0] to employ a Hubbard-Stratonovich transformation, in which the collective 
fields are introduced into the functional integral by inserting the following identities 



1 = Aq / D(p q exp 



d 4 p 



(3) 



N 2 J DA*DA q exp 

1 /A" 
G 2 V 2 



d 4 p 



A *Ct 
9 



G 2 



(4) 
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The quadratic terms in equations (3) and (4) cancel the original four-fermion interaction. 
As a result, the four-fermion interaction is cast into a Yukawa interaction between fermions 
and collective fields and a mass term for the collective fields. The fermionic fields then 
appear only quadratically and can be integrated out exactly, leaving an effective action for 
the bosonic collective fields alone. Having done the Gaussian integral for the fermions, the 
thermodynamic potential Q q reads 

^[0 ? ,A 9 ;T,//] = + ±-A 2 q + B> 

4lndetC 7 ^ + 7 >- ( ^ + ^ r .„_ *\ 9q)V 

(5) 



2 V K [7% + 7V0 - K + (j) q )\ 



r 



Disregarding the 7 matrix, the matrix of the last term of the above equation is 12 x 12. A 



q 



is a 6 x 6 matrix whose matrix elements are (A q ) ! °' = Cj 5 A q EijE 3a P , where i, j are flavor 
indices and a, j3 are color indices. 

The thermodynamic potential can be obtained by using imaginary time Green function 
method, with p = (2n + 1)^/(3 (/3=1/T). After performing the sums of the Matsubara 
frequences, the thermodynamic potential Q is written as 

T=u,d 

+Tln (1 + e -^ + ^) + + A 2 + 9 (C T _) y^ 2 - + A 2 

+2Tln (l + e^v 7 ^^) + 2Tln (l + e -l»&-)y/%~t&} } + B\ (6) 

where ( T ± = E T ± /i T , E T = a/ m* 2 + p 2 with m* = m T0 + <fi and 8(x) is a step function 
defined as 

The quark condensate (V>V) is related to via = —2Gi(^^)). 

In the four fermion interaction, form factors are not included [5|. Following Ref. 0], we 
introduce a cutoff, A g to regulate the above integral. After minimizing the thermodynamic 
potential with respect to and A, the following gap equations can be obtained. 

ml = m T0 + AG, [ { 1 - n F (E T ) - n F (E T ) + Cr+ [1 - 2n d (C T+ )] 



_ p 2 ^ j A T _ g (Cr-)A r _ ^ , 



7T" 
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where np(E T ) and uf(E t ) are the quark and anti-quark distributions, respectively, expressed 
as 

n F {E T ) = {exp [{E T - p T ) /T\ + l}" 1 , (10) 
n F (E T ) = {exp i(E T + p T ) /T] + l}" 1 . (11) 
The function rid(( T +) and nd(Cr-) are defined as 



nd((r-\ 



exp 



n d (Cr- 



exp 



0((r-)J( 2 T - + Al/T 



-1 



1 



-1 



(12) 



(13) 



From the thermodynamic potential Q, the pressure and quark density can be obtained as 
Pq = —Q and p T = dVt/dp T . One can define baryon density as p^ = | (p u + Pd)- In the 
chiral limit, the current quark mass m U Q=mdo — 0. In this case, equation (8) has a trivial 
solution m* = 0. The gap equation for A q also has a trivial solution A = 0. Therefore, for 
the massless quark, at some range of chemical potential or baryon density, there exist four 
sets of solutions: A. m* ^ 0, A ^ 0; B. m* ^ 0, A = 0; C. m* = 0, A / 0; D. m* = 0, 
A = 0. The system favors being in the phase with lowest thermodynamic potential. By 
comparing the thermodynamical potentials, one finds that the system will be in the phase 
with m* ^ 0, A = at low chemical potential (temperature) and in the phase with m* = at 
large chemical potential (temperature). However, at low chemical potential (temperature), 
the system cannot be in the quark phase but in the hadronic phase. Therefore, we need to 
replace the EOS for chiral symmetry breaking quark matter by that for hadronic matter. 



III. NUCLEAR MATTER 



For hadronic matter it is possible to work from the quark level, as in the quark meson 
coupling model or even in the NJL model. However, in the NJL model it is essential to 
incorporate the effect of confinement (for example through proper time regularization) in 
order to avoid a chiral collapse [15]. In the present work we adopt the more common hybrid 
approach and use the SU(3) quark mean field model in the hadronic phase. In that model 
the total effective Lagrangian is written: 



C 



(.•(I' 



£q0 + £>qM + £s£ + &VV + ^xSB + ^Am, + £-h, + £ 



(14) 



where C q o = q i^^du q is the free part for massless quarks. The quark-meson interaction £ 



qM 



-qO 

can be written in a chiral SU (3) invariant way as 

C qM = g s {^lM^ r + f fi M + f L ) - g v {^lI^I^l + *K7%*fl) 
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In the mean field approximation, the chiral-invariant scalar meson £ss and vector meson 
Cyv self-interaction terms are written as 

£ ES = -\ hx 2 {° 2 + C 2 ) + h (a 2 + c 2 ) 2 + h K + (A + h X °\ 



hx A ~ Vln4 + ^X 4 ln^, (16) 
4 Xo 3 ^oCo 



£vv = ~ — 2 (mluj 2 + m 2 p p 2 + m 2 ^ 2 ) + g A (u A + Qlu 2 P 2 + p 4 + 20 4 ) , (17) 
1 Xo 



r - x 



where 5 = 6/33; ao, Co an d xo are the vacuum expectation values of the corresponding mean 
fields a, ( and x- 

The Lagrangian £ x sb generates nonvanishing masses for the pseudoscalar mesons 

" ImlF^a + (V2m 2 K F K - ^F v ) A , (18) 
Xo L V V2 / . 

leading to a nonvanishing divergence of the axial currents which in turn satisfy the partial 
conserved axial-vector current (PCAC) relations for tc and K mesons. Pseudoscalar, scalar 
mesons and also the dilaton field x obtain mass terms by spontaneous breaking of chiral 
symmetry in the Lagrangian (1161) . The masses of u, d and s quarks are generated by the 
vacuum expectation values of the two scalar mesons a and 

In the quark mean field model, quarks are confined in baryons by the Lagrangian C c = 
— ^Xc^- Xc{ r ) is a confinement potential, i.e. a static potential providing the confinement 
of quarks by meson mean-field configurations. In the numerical calculations, we choose 
Xc( r ) = \k c r 2 , where k c = 1 GeV fm -2 . The Dirac equation for a quark field ^ under the 
additional influence of the meson mean fields is given by 

-ia ■ V + (3xc(r) + ftm* % = e*%, (19) 

where a = 7°7, ft = 7 , the subscripts i and j denote the quark i {i = u, d, s) in a baryon 
of type j (j = N, A, E, H). The quark mass m* and energy e* are defined as 

m* = -gia - g[C + m t0 (20) 

and 

e* = e t - glu - gi<P , (21) 

where is the energy of the quark under the influence of the meson mean fields. The 
effective baryon mass can be written as 

M* = Y,^<~ E l (22) 

i 

where n%j is the number of quarks with flavor "i" in a baryon j and E® is adjusted to give 
a best fit to the free baryon mass. 

Based on the previously defined quark mean field model the Lagrangian density for nuclear 
matter is written as 

C = $<yf% - M* N )^ + l -d^a + i^C^C + \d»xd»X - - X - Pilvf f» 

-guiny&u* - 9 P i'Bi^r 3 ^p fi + c M , (23) 
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where 

= d p u u - dyUJf, and p pu = d^p v - d u p p . (24) 

The term Cm represents the interaction between mesons which includes the scalar meson 
self-interaction the vector meson self-interaction Cyy and the explicit chiral symmetry 
breaking term C x sb, all defined previously The Lagrangian includes the scalar mesons a, £ 
and Xi an d the vector mesons uj and p. The interactions between quarks and scalar mesons 
result in the effective nucleon mass Mjy, The interactions between quarks and vector mesons 
generate the nucleon- vector meson interaction terms of equation (|23|) . The corresponding 

vector coupling constants gr w and g p are baryon dependent and satisfy the SU (3) relationship: 
Q v = -n"- — l a p = I i 

At finite temperature and density, the thermodynamic potential is defined as 

n = -h*L y f°° ^ f ln ^ + e - ( E* N{ k)-u N) /k B T^ + ln (! + e -(E h{k)+ u N ),k B Tx \ _ £ , ?5) 



where E* N {k) = y Mff + k 2 . The quantity z//v is related to the usual chemical potential, 
p N , by u N = p N — g^u — p. The energy per unit volume and the pressure of the system 

are respectively e = Cl — + vnPn and pu = — fi, where p^r is the baryon density. 

The mean field equation for meson 0j is obtained by the formula dVt/dipi = 0. For 
example, the equations for a, ( are deduced as: 

25 y 2 

k oX 2 a - Ak, (a 2 + C 2 ) a - 2k 2 a 3 - 2k 3X v( ~ ^X* + h m l F * 

XV 2 9m w 2 dm P, dM N,Ji\ n ^ 



A; oX 2 C - 4fci (a 2 + C 2 ) C - 4fc 2 C 3 - fc 3 va 2 - ^v 4 + |j ( V^m 2 ,^ - = 0, (27) 
where 

1 f°° k 2 M* 

W) = ^ / dfc ^y [nn(A:) + ™ n(A:) + np ^ k) + ™ p(A:)] • (28) 

In the above equation, n^{k) and n^(k) are the nucleon and antinucleon distributions, 
respectively, expressed as 

= UT?*(1\ 1 S // ^1 i 1 ( 29 ) 

exp [(£*(«) - v N ) /k B T\ + 1 

and 

exp l(E*(k) + u N ) /k B T\ + 1 

In the vacuum, the thermodynamic potential can have another solution where a = 0, ( = 
and x = 0- The energy difference of these two kinds of vacua is 

AE = C M {(J = 0, C = 0, X = 0) - C M {a = a , ( = Co, X = Xo)- (31) 
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It is interesting that the calculated value AE is close to the bag constant of quark matter. 
The numerical values of AE are listed in Table II. 

In the same way, the equations for the vector mesons, u and p, can be obtained: 

^mluj + Ag 4 u 3 + Ylg^p 2 = g"(p P + p n ), (32) 
Xo 

—m 2 p p + Ag A p 3 + I2g 4 u 2 p = -g^(p p - p n ), (33) 
Xo 6 

where p p and p n are the proton and neutron densities, expressed as 

1 f°° 

Pn = ^j dkk 2 [n q (k) - n q {k)) {N = p,n). (34) 
n Jo 

In nuclear matter, when the density is high, the volume of the nucleon is important. The 
nucleon cannot be treated as a point-like particle. The excluded volume effect has been 
considered in Ref . 25] . The key point is that the pressure of a system with excluded volume 
is the same as that of a point particle system, p B {T,p), if the chemical potentials of these 
two systems have the following relationship 

H = p + v p B (T,p), (35) 

where vq is the volume of a nucleon. pb(T,/i) is the pressure generated by baryons, i.e. 
Pb(T,£l) — Ph — Cm- The baryon density, entropy density and the energy density of the 
system (without meson contribution) are given by the usual thermodynamic expressions. 
We find the following relations: 

= ( ?PR\ - ( HER] ?E = P B (36) 
\ dp J T \ dp J T dp 1 + v op B ' 

c fdp B \ (dp B \ (dp B \ ( dp\ S B 

SB -\W)-\ln), + \^)Aw)-T^F B > (37) 

e B = TS B -p B + pp B = — — . (38) 

1 + v p B 

When the effect of the finite volume of the nucleon is included, the scalar density, (ipip), in 
equation (26) becomes 

1 f°° k 2 M* 

(W) = 2n 7 / dk—^- [n n {k) + n n (k) + n p (k) + n p {k)] . (39) 

n 2 (l + v p B ) J E*(k) 

The phase transition from hadronic matter to quark matter is determined by the following 
conditions: 

Ph=Pq, P P = 2fi u + p d , p n = p u + 2p d . (40) 

We should mention that at finite temperature, gluons make an important contribution to 
the pressure of quark matter. The gluon contribution to the pressure of quark matter is 
expressed as 

8 f°° k 3 

^=3Wn dk 7W^l- (41) 
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IV. NUMERICAL RESULTS 



In the NJL model, A q is chosen to obtain the quark condensate (qq) = —250 MeV. G\ is 
determined by the bold quark mass in vacuum m q = 350 MeV. There is some uncertainty 
in the coupling constant G 2 . The instanton interaction suggests GijG\ = 1/2, while in 
Ref. [5j, this ratio is 3/4. In our calculations we choose G^jG\ = 0.6, which is close to 
that of Ref. The parameter B' is determined by choosing the bag constant of quark 
matter to be 209MeV, i.e. pQ{m* — 0, fi q = 0) = — (209MeV) 4 . The numerical values of 
these three parameters are: A q = 634.7 MeV, G x =5.6 GeV~ 2 and (B 1 ) 1 ^ = 605.7 MeV. 
As for the nuclear matter phase, the parameters in the Lagrangian, which are shown in 
Table I, are determined by the saturation properties and the meson masses in vacuum. The 
corresponding nuclear matter properties and scalar meson masses are listed in Table II. 

In Fig. 1 we plot the nonzero solutions for the effective quark mass and energy gap 
versus baryon density at zero temperature. From the figure, one can see that the effective 
quark mass decreases continuously to zero with the increasing density. The energy gap also 
changes continuously with the density. However, at low temperature, the effective quark 
mass cannot decrease to zero continuously, as shown in Fig. 1. The energy gap cannot 
change continuously either - there is a first order phase transition. We have mentioned in 
section II that there are four kinds of solutions to Eqs. (8) and (9). At a given chemical 
potential the system will be in the phase with the lowest thermodynamic potential (highest 
pressure). In Fig. 2 we plot the pressure versus quark chemical potential for the four kinds 
of solutions at zero temperature. The solid, dashed, dotted and dash-dotted lines are for 
m* = 0, A = (phase D); m* = 0, A ^ (phase C); m* ^ 0, A ^ (phase A) and 
m* 7^ 0, A = (phase B), respectively. At small chemical potential, the system will be in 
the phase B: A = 0, m* ^ 0. At large chemical potential, the system will be in the phase 
C: A / 0, m* = 0. The phase transition takes place when the chemical potential is about 
342 MeV. Because this chemical potential is smaller than the effective quark mass in phase 
B, the density of this phase is zero. The effective mass and the energy gap versus chemical 
potential are shown in Fig. 3. At the transition point, fi q = 342 MeV, the changes in these 
quantities are discontinuous. 

With increasing temperature the energy gap will decrease. Fig. 4 shows the energy gap 
versus density at different temperatures. When the temperature is larger than about 52 MeV 
the Cooper pairs cannot form and the energy gap disappears. Fig. 5 shows the corresponding 
phase diagram. The solid and dashed lines indicate first and second order phase transitions, 
respectively. At low temperature, the transition from phase B to phase C is of first order. 
The transition from phase C to phase D is second order. This means that the two p — fi 
curves, corresponding to phase C and D, turn into a single curve at the transition point. 
We should mention that the result of the second order phase transition is an artifact of 
the mean field approximation. Early in 1974, Halperin, Lubensky and Ma have shown that 
the superconducting to normal phase transition can not be the second order 271 ] . In fact, 
this phase transition must be of the first order when gauge field fluctuation is included 
111, 3, 30, 31]. At high temperature the effective mass drops to zero continuously. The two 



p — fi curves of phase B and D turn into one curve smoothly and they do not cross each 
other as at low temperature. 

However, this phase diagram is for the transition from quark matter to quark matter, not 
for the transition from nuclear matter to quark matter. In nuclear matter the "elementary 
particles" are protons and neutrons. The equations of state for quark matter and nuclear 
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matter are different. The energy per baryon for nuclear matter is shown in Fig. 6, where 
nucleon volume effects are considered. The three curves correspond to the nucleon radius R 
= 0, 0.6 fm and 0.75 fm, respectively. The inclusion of finite volume effects makes the EOS 
for nuclear matter harder and generates a strong repulsive interaction between nucleons. 

The p — fi curves for nuclear matter and quark matter at zero temperature are shown 
in Fig. 7. The three solid curves are for nuclear matter and the two dashed lines for quark 
matter. If the nucleons are treated as point particles, there is no phase transition between 
nuclear matter and quark matter in this model. The chemical potentials for the nucleons are 
always smaller than the corresponding sums of quark potentials. When the volume effects 
are included, the nucleon chemical potentials will be increased and there is a first order phase 
transition between color superconducting quark matter and nuclear matter. The T — \x phase 
diagram for the transition between nuclear matter and quark matter is shown in Fig. 8. The 
solid and dashed lines indicate first and second order phase transitions, respectively. The 
letter H indicates the hadronic phase of nuclear matter. When the chemical potential is 
zero, the nuclear matter will turn into quark matter at temperature T = 148 MeV. At 
low temperature, the phase transition from nuclear matter to color superconducting quark 
matter will take place when the nucleon chemical potential, [in, is about 1.5 - 1.7 GeV. 

Figs. 9 and 10 show the T — p phase diagram for R = 0.6 fm and 0.75 fm, respectively. 
Since the phase transition is of first order, the densities of nuclear matter and quark matter 
are different. There is a coexistence state, M, in the figures. For the case of R = 0.6 fm, 
at zero temperature, when the nuclear density is about 0.76 fm -3 , quark matter starts to 
appear with density around 1.17 fm -3 . For the case R = 0.75 fm, the corresponding densities 
of nuclear matter and quark matter are 0.44 fm~ 3 and 0.85 fm~ 3 . From our calculation, we 
conclude that at low temperature, when the nuclear density is of order 2.5po - ^>Po, where 
Po is the saturation density of symmetric nuclear matter, there will be a phase transition. 
Since the phase transition between color superconducting quark matter and normal quark 
matter is of second order, the density remains the same during the phase transition. 

In our calculation, most parameters are fixed by well known physical quantities. Besides 
the nucleon radius R, there is only one parameter, the bag constant, which will sensitively 
affect the transition density and temperature. We discuss how the bag constant affects the 
transition temperature at zero baryon density. From the Fig. 11, one can see that with 
increasing temperature, the pressure of nuclear matter at zero baryon density increases very 
slowly. Until the temperature is around 180 MeV, the pressure remains low. However, for 
quark matter, the pressure increases very fast when the temperature is larger than about 
90 MeV. This is because the quark mass is much smaller than the nucleon mass. If the bag 
constant lies in the range 180 < B 1 ^ < 230 MeV, the transition temperature at zero baryon 
density is in the range 130 < T < 170 MeV. In other words, if the transition temperature is 
about 150 MeV the bag constant will be around (210 MeV) 4 which means that our choice 
of bag constant is reasonable. 

V. SUMMARY 

We have investigated the phase transition from nuclear matter to quark matter. The 
SU(3) quark mean field model, including nucleon volume effects, has been applied to get 
the EOS of nuclear matter. Color superconducting quark matter has also been studied, 
using the NJL model. The parameters in the models are determined by well known physical 
quantities. We now summarize the main results. 
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The phase transition between nuclear matter and quark matter is of first order. There is 
no critical temperature above which the phase transition is of second order. The so called 
second order phase transition is a transition from chiral symmetry breaking quark matter 
to chiral symmetry restored quark matter. However, in the real world, at low density or 
temperature, there is no chiral symmetry breaking quark matter, but rather one has nuclear 
matter. Therefore, the phase transition from nuclear matter to quark matter is always a 
first order transition, since the EOS of nuclear and quark matter are different. 

The effects of the finite volume of the nucleon are very important, especially at high 
density. The volume effects generate a strong repulsive interaction between the nucleons and 
increase their chemical potentials. The transition density is sensitive to the nucleon radius. 
In our numerical calculations, we chose R = 0.6 fm and 0.75 fm. At low temperature, when 
the nuclear density is of order 2.5p - §Po, the color superconducting quark phase will appear 
with a density 0.4 fm -3 larger than that of nuclear matter. 

Besides the nucleon radius, the phase transition is also sensitive to the bag constant B. 
If the bag constant is in the range 180 MeV < Z? 1 / 4 < 240 MeV, the transition temperature 
will be 130 MeV < T < 170 MeV. The bag constant B can be explained as the energy 
difference, AE, between two kinds of vacua. In the 57/(3) quark mean field model, the 
calculated AE is reasonable, giving a transition temperature around 150 MeV. 

At low temperature, when the density increases, the nuclear matter will always turn into 
color superconducting quark matter. When the temperature is larger than about 52 MeV, 
the superconductivity phase disappears. The transition from normal quark matter to color 
superconducting quark matter is of second order, while the transition from nuclear matter 
to quark matter is of first order. 

In conclusion, we note that our discussions concerning the order of the phase transitions 
investigated here have been based upon mean field theory. In fact, for the transition from 
superconducting quark matter to normal quark matter, with the gluon field fluctuation, this 
transition must be of the first order. It clearly will be important to explore such effects 
in the current model. In addition, we have presented results without any consideration of 
strangeness. As strange quarks may be expected to appear at perhaps 2.5 to 3 times normal 
nuclear matter density (for matter in beta equilibrium, as in a neutron star), it will also be 
of interest to extend the present model to include strange quarks. 
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TABLE I: Parameters of the SU (3) quark mean field model. 



R (fm) 


k 


h 


k 2 


^3 


k 4 


9s 


9v 


9a 





3.97 


2.18 


-10.16 


-4.15 


-0.14 


4.76 


8.70 


15.0 


0.6 


3.95 


1.83 


-10.16 


-4.56 


-0.10 


4.76 


8.70 


15.0 


0.75 


3.94 


1.40 


-10.16 


-5.13 


-0.03 


4.76 


8.70 


15.0 



TABLE II: Nuclear matter properties and scalar meson masses in the SU(3) quark mean field 
model. 



R (fm) 


Po (fm" 3 ) 


E/A (MeV) 


M* N /M N 


K (MeV) 


(AE 1 ) 1 / 4 (MeV) 


m a (MeV) 


m c (MeV) 





0.16 


-16.0 


0.742 


303 


220.2 


487.8 


1168.0 


0.6 


0.16 


-16.0 


0.731 


359 


215.5 


490.0 


1173.4 


0.75 


0.16 


-16.0 


0.713 


473 


209.5 


491.6 


1181.2 



< 200 

* 

E 



' 1 1 1 


i | ■ | i 


\m* 






A 





0.0 0.2 0.4 0.6 0.8 1.0 

P B (fm- 3 ) 

FIG. 1: Effective quark mass and superconducting energy gap calculated in the NJL model versus 
baryon density at zero temperature. 
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FIG. 2: Thermodynamic potential of quark matter versus quark chemical potential, fj, g , at zero 
temperature. The solid, dashed, dotted and dash-dotted lines are for m* = 0, A = 0, m* = 0, 
A / 0, m* / 0, A / and m* / 0, A = 0, respectively. 
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FIG. 3: Effective quark mass and superconducting energy gap versus quark chemical potential at 
zero temperature. 
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FIG. 5: T-fi phase diagram for quark matter. The solid line stands for the first order phase 
transition and dashed lines stand for the second order phase transition. B, C and D correspond to 
the different phases explained in the text. 
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FIG. 6: The energy per baryon of nuclear matter versus baryon density. The radius of the nucleon 
is chosen to be 0, 0.6 fm and 0.75 fm, respectively. 
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FIG. 7: Pressure of nuclear matter and quark matter versus chemical potential. The solid lines 
are for nuclear matter and dashed lines are for quark matter. 
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FIG. 8: T-fi phase diagram for nuclear matter and quark matter. The solid lines are for first order 
phase transition between nuclear matter and quark matter. The dashed line is for the second order 
phase transition between superconducting quark matter and normal quark matter. H is for the 
phase of nuclear matter and C and D are for the quark phases. 
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FIG. 9: T-p phase diagram for nuclear matter and quark matter. The radius of nucleon is chosen 
to be 0.6 fm. H is for the phase of nuclear matter. M is for the mixed phase. C and D are for the 
quark phases. 
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FIG. 10: T-p phase diagram for nuclear matter and quark matter. The radius of the nucleon is 
chosen to be 0.75 fm. H is for the phase of nuclear matter. M is for the mixed phase. C and D are 
for the quark phases. 
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FIG. 11: Pressure of nuclear and quark matter at zero baryon density versus temperature. The 
solid lines are for quark matter with the different bag constants. The dashed line is for nuclear 
matter. 
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